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We report tunnel diode resonator measurements of in-plane London penetration depth, A(T), 
in optimally-doped single crystals of Fe(Teo.5gSeo.42) with Tc ^ 14.8 K. Systematic measurements 
were carried out for six samples with different size and surface roughness. The power-law behavior, 
AA(r) = AT" was found for all samples with the average exponent Uavg ~ 2.3 ±0.1 and the 
pre-factor Aavg = 1.0 ± 0.2 nm/K^'^. The average superfluid density is well described by the self- 
consistent two-gap 7 model resulting in Aj{0) /ksTc = 1.93 and Ajj{0)/kBTc = 0.9. By analyzing 
the data obtained on samples of different size and deliberately introduced surface roughness, it is 
concluded that the calibration procedure used to obtain A(r) from the measured TDR frequency shift 
is quite robust and the uncertainty in sample dimensions and the nature of surface roughness play 
only minor role. The exponent n, directly related to the superconducting gap structure, remains 
virtually unchanged. The calibration - dependent pre-factor A shows some variation, but stays 
within reasonable margin ruling out some recent suggestions that surface conditions can significantly 
affect the results. Our results confirm that precision global measurements provide the most objective 
information regarding temperature - dependent London penetration depth. 



The crystals of iron-based superconductors are mica- 
ceous and soft, which prompts questions related to sam- 
ple structure and homogeneity. Many measurements as- 
sume perfect uniformity of the bulk properties and per- 
fectly flat surface. This is especially important for the 
measurements of London penetration depth, A(r), which 
is one of the key quantities that can be measured with 
great precision providing valuable information on the 
quasiparticle excitations, thus the superconducting gap 
structure. When the entire sample is measured, several 
factors may affect the outcome. Size and shape enter 
the cahbration factors that is used to convert raw data 
into the actual length. Chemical inhomogeneity may af- 
fect the measurements as an intrinsic source. One recent 
example of the experimental controversy is the compar- 
ison of A(r) estimated from the tunnel-diode resonator 
(TDR) and first critical field measurements in single crys- 
tals of Fe(Se,Te) [1]. In order to reconcile two data sets 
the authors had to assume that the effective sample di- 
mension, R, in the TDR calibration should be five times 
larger than the calculated from the sample dimensions, 
presumably due to surface roughness. 

In order to test the effects of sample size, shape 
and surface roughness we measured London penetration 
depth in optimally doped Fe(Teo.58Seo.42) with Tc ^ 14.8 
K using a TDR technique. A series of measurements were 
carried out on parts of the same sample with different 
sizes and surface roughness. We find a robust power-law 
variation, A(r) ^ AT" with n — 2.3 ± 0.1 among all mea- 
surements indicating intrinsic behavior. For the purpose 
of investigating how surface roughness affects AA(T), the 
edges of one of the samples were deliberately damaged by 
a blade, see Fig. [T] Measurements of this sample before 




FIG. 1. (Color online) (a) Sample 2 - before and after (b) 
(Sample 2-R) deliberate roughening, see Table ^ 



and after damaging the edges show no significant change 
of the pre-factor A or the exponent n. These results 
suggest that London penetration depth measured by a 
TDR technique represents an intrinsic quantity provided 
that the samples exhibit homogeneous superconductiv- 
ity. We emphasize that it is the exponent n that is most 
important quantity to characterize the gap structure and 
possible pair-breaking effects. 

The iron chalcogenide family of superconductors FeCh 
(Ch = Se/Te/S) was discovered in 2008 shortly af- 
ter the discovery of the first iron based superconductor 
La(Oi_a;Fa;)FeAs [3]. The structure of this family is sim- 
pler than those of other Fe-based superconductors (FeSC) 
[3H7] . The square planar sheets of Fe ions are alternated 
by the distorted tetrahedra of Chalcogen ions. Even with 
simpler structure, the electronic structure of FeCh is sim- 
ilar to other FeSCs. Superconductivity in FeCh appears 
upon substitution of Te with Se or S [SI |9], or control 
of the amount of excess Fe [21 E] • The superconducting 
transition temperature T^ is lower than those in most 



TABLE I. List of samples and their physical constraints. R 
is the effective dimension calculated from the dimensions of a 
sample [TOJ . 



Sample 


Dimensions (/im ) 


R (pm) 


Edge condition 


1 


712 X 488 X 40 


33.5 


clean cut 


1-A 


491 X 270 X 40 


22.6 


clean cut 


1-B 


489 X 443 X 40 


28.4 


clean cut 


2 


702 X 455 X 40 


32.8 


clean cut 


2-R 


680 X 430 X 40 


31.5 


rough edge 


2-R-C 


527 X 421 X 40 


28.6 


clean cut 



other FeSCs. In Fey{Tei-xSex) at optimal doping, Tc 
reaches a value of ~ 14 K and '--^ 36 K under high pres- 
sure [H [T2HT5] . A series of theoretical and experimental 
studies suggest that superconductivity in FeCh could be 
magnetically mediated following the spin-fluctuation pic- 
ture [I7H2Q]. 

The pairing mechanism of FeC/i has been explored 
from various measurements [211 I22j . In nuclear mag- 
netic resonance (NMR) study of polycrystalline FeSe {Tc 
= 8 K), the unconventional pairing was suggested from 
the absence of the NMR coherence peak in conjunction 
with the power-law temperature dependence of the spin- 
relaxation rate, 1/Ti ~ T^ [5T]. This can be explained 
with the unconventional order parameter such as a fully 
gapped s± or a nodal gap pairing. In thermal conductiv- 
ity measurements on FeSe^; single crystals, the multi-gap 
nodeless superconductivity was suggested [22] , 

The London penetration depth in FeCh has been mea- 
sured by using different methods pQ E51 - E5] . The result 
from muon spin rotation (/^SR) in FeSe^; is consistent 
with either anisotropic s-wave or a two-gap extended s- 
wave pairing (23j . Radio- frequency TDR measurements 
of Aa6(T) in Fey(Tei_^Se^) (x = 0.37 ^, a; = 0.44 [15], 
and X =0.45 [1 ) by three different groups found simi- 
lar power-law exponent n and some variation in the pre- 
factor A: n ~ 2.0, A -3.7 nm/K" [Ml, n ~ 2.2, A ~ 
0.9 nm/K" [25] and n - 2.0, A - 4 nm/K" [1]. The 
full-range temperature behavior was also very similar. 
However, as mentioned above, an attempt to compare 
the TDR results with A(r) estimated from the measured 
first critical field, Hd using Ginzburg-Landau formula 
(valid only at Tc) resulted in a large discrepancy forcing 
the Authors to suggest a very different effective dimen- 
sion of i? = 70 /im as compared to the one calculated 
from the sample dimensions, i? = 14 /um Jj. 

Single crystals of Fe(Tei-.a;Se2;) were synthesized us- 
ing a flux method as reported before [13]. The samples 
studied in this work had composition of Fe : Te : Se 
= 1.00 ± 0.01 : 0.58 ± 0.01 : 0.42 ± 0.01 determined 
from the energy dispersive x-ray spectroscopy (EDXS). 
Initially two samples 1 and 2 were prepared by careful 
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FIG. 2. (Color online) Temperature dependent frequency 
shift, A/(r) for each sample listed in Table [l] Tc is consis- 
tently determined to be 14.8 K by using onset curves at the 
phase transition. Inset : A/(T) normalized based on the value 
at 16 K. The high temperature region is zoomed in showing 
the phase transition. 



cleaving and cutting processes. Then Sample 1 was split 
into Samples 1-A and 1-B. To investigate the influence 
of the edge roughness, deliberate damage was done to 
Sample 2 by a razor blade. The original (Sample 2) and 
roughened ones (Sample 2-R) are shown in Fig.fl] To pro- 
duce heaviest roughness possible, the damaging process 
was repeated multiple times resulting in a loss of about 
10 % of volume. After the measurements, the rough sur- 
face of Sample 2-R was cut again as clean as possible to 
remove the roughness and the cleaned one was labeled 
as 2-R-C. The penetration depths was measured at each 
stage of the described procedures. The dimensions and 
edge condition of samples are summarized in Table [l] 

The in-plane London penetration depth, A(T), was 
measured using a self-oscillating tunnel-diode resonator 
(TDR) technique [TUl IMl [27] . A sample under study is 
mounted on a sapphire rod and inserted into an inductor 
coil of a LC tank circuit. To measure the in-plane pene- 
tration depth, the sample is placed with its c-axis along 
the direction of ac-field (Hac) induced by the inductor 
coil. Since Hac ^ 20 mOe is weak enough (^ i?ci), the 
sample is in the Meissner state, so its magnetic response 
is determined by the London penetration depth. The fre- 
quency shift in TDR, A/ = /(T) ~ /o, is used to obtain 
the magnetic susceptibility x(2^) from A/ = — G47rx(T). 
Here /o = \j2'n\fLC ^ 14 MHz is the resonant frequency 
of an empty resonator, G = /oT4 7214(1 — A^) is a geo- 
metric factor defined by the coil and sample volumes, Vc 
and Vs , and N is the demagnetization factor. Calibration 
constant G is directly measured by pulling the sample out 
of the coil at the lowest temperature. In the Meissner 




FIG. 3. (Color online) Temperature dependent penetration 
depth, AA(T) for each sample, (a) High temperature region 
and (b) Low temperature region of AA(T). The arrows indi- 
cate the upper-limits of fitting. The dashed lines are repre- 
sentative fits for each sample conducted up to Tc/?). All the 
fitting results are summarized in Fig. |4] 
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FIG. 4. (Color online) Results of the power-law fits to AA(T) 
= AT". Four different upper-limits, indoicated by arrows in 
Fig. [slb), were used. Panel (a): the exponent n, obtained by 
keeping A and n as free parameters. Panel (b): the pre-factor 
A, obtained at a fixed Uavg ~ 2.3, which is the average among 
all seven samples shown in Panel (a). 



state, A can be obtained from the magnetic susceptibil- 
ity, Xi following the relation Anx = (A/i?)tanh(i?/A) — 1 
where R is the effective dimension of the sample [TU] . 

Figure [2] shows the raw data, temperature dependent 
frequency shift, A/(T), for samples listed in Tabic |lj It 
is obvious that all the samples have quite sharp and well 
overlaping transition curves with the width of the super- 
conducting transition of about ^ 0.3 K, so it seems that 
intrinsic superconducting properties were consistently re- 
producible between the samples (all samples come from 
the same large piece). 

Above Tc the signal, A/(r), saturates due to two pos- 
sible reasons. One, so-called sample-size-limited regime, 
is due to the size of a sample, which affects the tanhi?/A 
factor in the expression for the differential susceptibil- 
ity. Another reason for the saturation, so-called skin- 
depth hmited, occurs if the penetration depth becomes 
limited by the normal state skin depth. At Tc, the esti- 
mate value of the skin depth is about twice of that es- 
timated from our calibration for X{T) 28 . Thus, the 
samples under current study is in sample-size-limited 
regime, so their penetration depth shifts, AA(r), are cal- 



culated from A/(T) in Fig. [s] following the calibration 
procedure explained above. Figure pfa) shows full tem- 
perature range which reveals sharp phase transitions at 
14.8 K. The low temperature behavior of AA is shown 
in Fig. Islb). From these calibrated values we see that 
the penetration length in the normal state is compara- 
ble to the values of R, so we are in a sample-size limited 
regime. Therefore, the total frequency shift between the 
lowest temperature (~ 500 mK) and the normal state 
temperature (T > T^), A f total, is related to the sample 
dimension. For instance, the largest Sample I has the 
biggest Af total and the smallest Sample 1-A the smallest 
A f total- The inset shows Af(T) normalized by the value 
at 16 K. 

For quantitative analysis of the low temperature be- 
havior, a power-law fit, AA(r) = AT" was performed for 
all six samples. To check how robust is the power law, in 
each case, four different upper temperature limits for the 
fit were used T^/S, Tc/4:, T^/b and T^/Q, shown schemati- 
cally by arrows in Fig.l3|b). The results of the fitting are 
summarized in Fig. [4] Clearly, the fit coefficients remain 



fairly constant with small (and expected) deviations at 
the lowest and the highest limits. Below we discuss the 
results for Tc/A chosen as the upper limit for the fit. 

As shown in Fig. H (a), in general, the exponent n is 
rather steady among all cases even though there are small 
variations. From the best fits with T /Tc upper limit, the 
average exponent Uavg — 2.3 ±0.1 is quite comparable 
to the previous reports of n =2.1 |21j, 2.2 [351 and 2.0 
[1]. In conjunction with the superfluid density analysis 
discussed later, the exponent n = 2.3 can be explained 
by the nodeless two-gap pairing symmetry with strong 
pair breaking effect. The behavior of the pre- factor A is 
summarized in Fig. B[b) . Clearly, there is a very weak 
dependence on the fitting range, but there is a more sub- 
stantial change between the samples. The average over 
all samples pre-factor Aavg = 1.0±0.2 nm/K^'^ is similar 
to other reports [MIUS]- 

Let us now try to understand the observed variation of 
A for different samples. We choose the values obtained 
from the fitting range up to Tc/4 where the fit quality 
was the best for all samples. We started by measuring 
Sample 1 that showed A — 0.9 nm/K^'^. After then 
Sample 1 was cut into Samples 1-A and 1-B. Every at- 
tempt was made to cut as clean as possible, but the pre- 
factor A has increased for both samples 1-A and 1-B to 
1.3 nm/K^'^. Although this change is not too significant, 
this result suggests that the effect of the rough edges be- 
comes more pronounced as the sample size decreases. To 
estimate the effect of roughness itself when the size of a 
sample remains the same, Sample 2-R was made out of 
Sample 2 by deliberately damaging the edge with a ra- 
zor blade. However, it turns out that the increment of A 
was only by 0.4 nm/K^-^ from 0.7 nm/K^-^ for Sample 2 
to 1.1 nm/K^-^ for Sample 2-R. This increment is simi- 
lar to that between Sample 1 and two Samples 1-A and 
1-B. This means that the intentional roughness created 
by coarse roughening a sample doesn't significantly affect 
the result and it also means that the effective sample di- 
mension R used to calibrate the TDK data (so that it 
directly affects the pre-factor A) is very close to the cal- 
culated value. Otherwise, coarse roughening would make 
it many times larger due to the increased surface area ex- 
posed to the field. Furthermore, we tried to remove the 
rough surface from Sample 2-R by cutting the edges as 
clean as possible. The cleaned Sample 2-R-C showed that 
A has decreased by 0.1 nm/T^-^ from that of Sample 2- 
R. This re-enforces our conclusion that surface roughness 
is not a dominant source determining A and plays only 
minor role. We also conclude that the explanation pro- 
posed to reconcile TDR and Hd measurements [T] does 
not hold the ground. Instead, we think that the field of 
first penetration measured by the Hall-probe was deter- 
mined by the surface barrier rather than if ci . Finally, we 
conclude that global measurements provide a reliable es- 
timate of the overall behavior of the London penetration 
depth by sampling all sample surfaces in contrast to the 
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FIG. 5. (Color online) Superfluid density ps{T) = 

(A(0)/A(r))^ versus T/Tc. A(0) = 560 nm obtained from the 
previous report [M]. (a) Ps{T) for all six samples, (b) The 
average superfluid ndensity for all six samples, pj"* (Red cir- 
cle), is fltted to the two-gap 7 model [52|. Black solid line is 
the fltting result with p = 7Pj- + (1 — j)Pij ■ The inset shows 
two gaps acquired from the fit, A^ and A^^. 



local probes that are affected by the surface topography 
on the mesoscopic scale. 

The superfluid densities ps — (A(0)/A(T))^ for all sam- 
ples are shown in Fig. [51 The absolute penetration depth, 
A(0) = 560 nm was determined by TDR measurements 
of Al-coated sample [21] ■ The variation of ps among dif- 
ferent samples is shown in Fig. p^ (a). For the purpose 
of fitting, the average p™^ is calculated from ps of sam- 
ples 1 and 2 which have lowest A values. The fitting 
was done with a self-consistent two-gap 7-model, where 
two gaps are calculated self-consistently at each temper- 
ature and at each iteration. The total superfluid density 
is given by p = jp^ + (1 — 'y)pjj ^9\. The partial den- 
sities of states are chosen to be equal on the two bands, 
ni ~ n2 = 0.5 and Debye temperature of 230 K was used 
to calculate the experimentally observed Tc = 14.8 K, 
which fixes the coupling constants (we used An). Fig- 



urelsFb) demonstrates a good agreement between experi- 
mental p™^ (symbols) and fitting (black solid line). The 
parameters acquired from the fit are: An — 0.66, A22 — 
0.44, A12 = 0.07, Ae// = 0.34, 7 = 0.75 and T, = 14.95 
K. This result indicates that 75 % contribution of super- 
fluid density comes from the band with p^ which has the 
larger gap A^-. The inset in Fig. [5] (b) shows the behav- 
iors of two energy gaps A^ and A^j- versus temperature. 
Clearly, the smaller gap has significantly non-BCS tem- 
perature dependence. The zero temperature energy gaps 
A,(0) and A„(0) are 2.5 meV (A,(0)/fcBTc = 1.93) and 
1.1 meV (A„(0)/fcBTc = 0.9) respectively. The super- 
conducting energy gaps has also been probed in different 
experiments. From /nSR ^301 |3T] and penetration depth 
[21] measurements, two s-wave energy gaps were reported 
with gap values similar to our results. /xSR studies in 
FeSeo.sTeo.s [301 ISI] revealed two gaps of A/arge ~ 2.6 
meV and A^man ~ 0.5-0.87 meV, and the penetration 
depth study [2^ also showed that Ajarge ~ 2.1 meV and 
As ~ 1.2 meV. According to scanning tunneling spec- 
troscopy study, only one s-wave gap A ^ 2.3 meV was 
observed in FeSeo.4Teo.6 [32 , which is similar to the large 
gap A/(0) of our result. However, rather large single or 
multi-gaps were reported from specific heat [33], opti- 
cal conductivity [34], point-contact Andreeve reflectiv- 
ity ,35,, and angle-resolved photoemission spectroscopy 
[55] suggesting strong-coupling superconductivity. The 
electronic specific heat in Fe(Teo. 57860.43) [33] revealed 
two energy gaps with A;arge '^ 7.4 meV and l^smaii ~ 
5.0 meV. From the optical conductivity in FeTco. 55860. 45, 
two large energy gaps were also found with A/^rge ~ 5.1 
meV and ^smaii ~ 2.5 meV. The point-contact Andreev 
reflectivity in FeTco. 55869. 45 is consistent with single gap 
s-wave symmetry with A (at 1.70 K) ~ 3.8 meV. Angle- 
resolved photoemission spectroscopy in FeTeo.78eo.3 [36] . 
an s-wave single gap of A ^ 4 meV was also observed. 
Overall, the pairing symmetry in FeTei_2:8ea; is still un- 
der debate. 



To summarize, London penetration depth was 
measured in optimally-doped single-crystals of 
Fe(Teo. 58 860.42) of different size, shape and surface 
conditions to separate extrinsic and intrinsic effects 
determining the temperature variation of A(T). Even 
though there is a natural dispersion among six different 
samples, the average exponent Uavg and pre-factor A^vg 
are found to be 2.3 ± 0.1 and 1.0 ± 0.1 nm/K^-^. The 
two-band 7 model fits the superfluid density rather well 
resulting in A,(0)//cbTc = 1.93 and A„(0)/fci3Tc = 
0.9. From comparison of six samples, it is found that 
there exists unavoidable micro(meso)-scopic surface 
roughness on the sample edges and the effect of the 



roughness on pre-factor A becomes more important 
as the sample size decreases, whereas the exponent n 
remains practically unchanged. Overall, we conclude 
that global measurements provide an accurate and 
objective way to determine London penetration depth 
in pnictide superconductors. 
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